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1 Introduction 



Arrays of weakly coupled Josephson junctions provide an experimental realiza- 
tion of the two dimensional (2D) XY model physics. A Josephson junction 
ladder (JJL) is the simplest quasi-one dimensional version of an array in a mag- 
netic field recently such a system has been the subject of many investigations 
because of its possibility to display different transitions as a function of the mag- 
netic field, temperature, disorder, quantum fluctuations and dissipation. In this 
contribution we address the phenomenon of fractionalization of the flux quan- 
tum in a fully frustrated JJL, the basic question being how the phenomenon 
of Cooper pair condensation can cope with properties of charge (flux) fraction- 
alization, typical of a low dimensional system with a discrete (Z 2 in our case) 
symmetry. Then we discuss how it is deeply related to the issue of topological 
order. 

We must recall that charge fractionalization has been successfully hypothe- 
sized by R. Laughlin to describe the ground state of a strongly correlated 2D 
electron system, a quantum Hall fluid, at fractional fillings v = 2 p+i > P = L 2, 
In such a system charged excitations are present with fractional charge (anyons) 
and elementary flux — . Furthermore the phenomenon of fractionalization of the 
elementary flux has been found more recently in the description of a quantum 
Hall fluid at non standard fillings v = m ™ +2 El El- within the context of 2D 
Conformal Field Theories (CFT) with a Z m twist. 

In Refs. [5| it has been shown that the presence of a Z 2 symmetry accounts 
for more general boundary conditions for the propagating electron fields which 
arise in quantum Hall systems in the presence of impurities or defects. Further- 
more such a symmetry is present also in the fully frustrated XY (FFXY) model 
or equivalently, see Ref. [E]PP, in two dimensional Josephson junction arrays 
( J JA) with half flux quantum | ^| threading each square cell and accounts for 
the degeneracy of the ground state. 

It is interesting to notice that it is possible to generate non trivial topologies, 
i.e. the torus, in the context of a CFT approach. That allows in our case to 
show how non trivial global properties of the ground state wave function emerge 
and how closely they are related to the presence of half flux quanta, which can 
be viewed also as "topological defects" . 

The concept of topological order was first introduced to describe the ground 
state of a quantum Hall fluid 7 . Although todays interest in topological order 
mainly derives from the quest for exotic non-Fermi liquid states relevant for 
high T c superconductors J5J, such a concept is of much more general interest [Jj]. 

Two features of topological order are very striking: fractionally charged 
quasiparticles and a ground state degeneracy depending on the topology of the 
underlying manifold, which is lifted by quasiparticles tunnelling processes. For 
Laughlin fractional quantum Hall (FQH) states both these properties are well 
understood jTH], but for superconducting devices the situation is less clear. 

Josephson junctions networks appear to be good candidates for exhibiting 
topological order, as recently evidenced in Refs. C2 by means of Chern- 
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Simons gauge field theory. Such a property may allow for their use as "pro- 
tected" qubits for quantum computation. 

The aim of this contribution is to show that the twisted model (TM) well 
adapts to describe the phenomenology of fully frustrated J JL with a topological 
defect and to analyze the implications of "closed" geometries on the ground state 
global properties. In particular we shall show that fully frustrated Josephson 
junction ladders (JJL) with non trivial geometry may support topological order, 
making use of conformal field theory techniques 1 . A simple experimental test 
of our predictions will be also proposed. 

The paper is organized as follows. 

In Section 2 we introduce the fully frustrated quantum Josephson junctions 
ladder (JJL) focusing on non trivial boundary conditions. 

In Section 3 we recall some aspects of the m-reduction procedure m 
particular we show how the m — 2, p — case of our twisted model (TM) P] 
well accounts for the symmetries of the model under study. In such a framework 
we give the whole primary fields content of the theory on the plane and exhibit 
the ground state wave function. 

In Section 4 the symmetry properties of the ground state conformal blocks 
are analyzed and its relation with their topological properties shown. 

In Section 5, starting from our CFT results, we show that the ground state is 
degenerate, the different states being accessible by adiabatic flux change tech- 
niques. Such a degeneracy is shown to be strictly related to the presence in 
the spectrum of quasiparticles with non abelian statistics and can be lifted non 
perturbatively through vortices tunneling. 

In Section 6 some comments and outlooks are given. 

In the Appendix we recall briefly the boundary states introduced in Ref. [H] 
in the framework of our TM. 

2 Josephson junctions ladders with Mobius bound- 
ary conditions 

In this Section we briefly describe the system we will study in the following, 
that is a closed ladder of Josephson junctions (see Fig.l) with Mobius bound- 
ary conditions. With each site i we associate a phase ipi and a charge 2erij, 
representing a superconducting grain coupled to its neighbours by Josephson 
couplings; rij and (pi are conjugate variables satisfying the usual phase-number 
commutation relation. The system is described by the quantum phase model 
(QPM) Hamiltonian: 

H = -^J2 (£) - E cos - n - > C 1 ) 

» (ij) 

where Ec — ^q— (C being the grain capacitance) is the charging energy at 
site i, while the second term is the Josephson coupling energy between sites i 
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Figure 1: Josephson junction ladder with Mobius boundary conditions 



and j and the sum is over nearest neighbours. Ay = |p f? A-dl is the line 
integral of the vector potential associated to an external magnetic field B and 
$o = is the superconducting flux quantum. The gauge invariant sum around 
a plaquette is J2 p Ay = 2nf with / = , where $ is the flux threading each 
plaquette of the ladder. 

Let us label the phase fields on the two legs with tp\ , a = 1, 2 and assume 
Eij = E x for horizontal links and E^ — E y for vertical ones. Let us also make 



the gauge choice Ay = +nf for the upper links, A 



-irf for the lower 



ones and Ay = for the vertical ones, which corresponds to a vector potential 
parallel to the ladder and taking opposite values on upper and lower branches. 
Thus the effective quantum Hamiltonian can be written as 0: 



« = 



d 



E 



d 



dip 



(2) 



E x £ cos (VS - + {-If nf) + E y cos - ^ (2) )( 



a=l,2 



The correspondence between the effective quantum Hamiltonian (J2J) and our 
TM model can be best traced performing the change of variables [2]: ip^ — 

(2) 

Xi + (pi, ipl = Xi - (pi, so getting: 



» = -¥e 



d 
dX* 



cos {X l+1 - X t ) cos - 4>i - nf) + E y cos (2&)] (3) 



where Xj, (i.e. , "ft ) are om Y phase deviations of each order parameter 
from the commensurate phase and should not be identified with the phases of 
the superconducting grains [2|. 
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When f = \ and Eq — (classical limit) the ground state of the ID 
frustrated quantum XY (FQXY) model displays - in addition to the continuous 
U(l) symmetry of the phase variables - a discrete Zi symmetry associated with 
an antiferromagnetic pattern of plaquette chiralities \y = il> measuring the 
two opposite directions of the supercurrent circulating in each plaquette. The 
evidence for a chiral phase in Josephson junction ladders has been investigated 
in Ref. ^3] while a field theoretical description of chiral order is developed in 



Performing the continuum limit of the Hamiltonian lpj|): 



H 



Ec 
2 



dx 



dx 



E,. 



(~ 

\dX 

dX 
dx 



d_ 



E y cos (24>) 



(4) 



we see that the X and cf> fields are decoupled. In fact the X term of the above 
Hamiltonian is that of a free quantum field theory while the (f> one coincides with 
the quantum sine-Gordon model. Through an imaginary-time path-integral 
formulation of such a model |16j it can be shown that the ID quantum problem 
maps into a 2D classical statistical mechanics system, the 2D fully frustrated XY 

model, where the parameter a — (j^j plays the role of an inverse temperature 
For small Ec there is a gap for creation of kinks in the antiferromagnetic 
pattern of \v an< ^ the ground state has quasi long range chiral order. We work 
in the regime E x 3> E y where the ladder is well described by a CFT with central 
charge c = 2. 

We are now ready to introduce the modified ladder £Q, see Fig. 1. In order 
to do so let us first require the <f^ a \ a — 1,2, variables to recover the angular 
nature by compactification of both the up and down fields. In such a way the 
XY-vortices, causing the Kosterlitz-Thouless transition, are recovered. As a 
second step let us introduce at point x — a, defect which couples the up and 
down edges through its interaction with the two legs, that is let us close the 
ladder and impose Mobius boundary conditions. In the limit of strong coupling 
such an interaction gives rise to non trivial boundary conditions for the fields 
5 . In the following we give further details on such an issue, in particular we 
adopt the m-reduction technique , which accounts for non trivial boundary 
conditions for the Josephson ladder in the presence of a defect line. In the 



Appendix the relevant chiral fields tp, 



(a) 



1,2, which emerge from such 



conditions, are explicitly constructed, by using the folding procedure. 



3 m-reduction technique 

In this Section we focus on the m-reduction technique for the special m = 2 case 
and apply it to the system described by the Hamiltonian (0J. In the Appendix 
each phase field tp^ is written as a sum of two fields of opposite chirality 
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defined on an half-line, because of the presence of a defect at x — 0. Within 
a " bosonization" framework it is shown there how it is possible to reduce to a 
problem with two chiral fields <pi a \ a = 1,2, each defined on the whole a;— axis, 
and the corresponding dual fields. Now we identify in the continuum such chiral 
phase fields <+pf\ a = 1,2, each defined on the corresponding leg, with the two 
chiral fields Q {a \ a = 1, 2 of our CFT, the TM, with central charge c = 2. 

In order to construct such fields we start from a CFT with c = 1 described 
in terms of a scalar chiral field Q compactificd on a circle with radius R 2 = 2, 
explicitly given by: 

Q(z) = q-iplnz + i 1 %2 — z~ n (5) 

with a n , q and p satisfying the commutation relations [a„,a„'] = n8 n _ n ' and 
[QtP] = i\ its primary fields are the vertex operators U a ' (z) =: e 1 ™ 1 ®^' : 
and cti = -^=, I = 1,2. It is possible to give a plasma description through 

the relation \^\ 2 = e-< m -st where ijj (zi, z N ) = (Na\ U a < (zi)\o\ = 

N — 

llj<j=i ( z i ~ z j) 2 ls the ground state wave function. It can be immediately 

seen that H e ff = —I 2 J2i<j=i m \ z i ~ z j\ an( i P = 'jjz = 1) tnat i s om y vorticity 
v = 1, 2 vortices are present in the plasma. 

Starting from such a CFT mother theory one can use the m-reduction proce- 
dure, which consists in considering the subalgcbra generated only by the modes 
in eq. (0) which are a multiple of an integer to, so getting a c = to orbifold CFT 
(daughter theory i.e. the twisted model (TM)) 0. With respect to the special 
m = 2 case, the fields in the mother CFT can be organized into components 
which have well defined transformation properties under the discrete Zi (twist) 
group, which is a symmetry of the TM. By using the mapping z — » z 1 / 2 and by 
making the identifications a 2n +i — > V2a n +i/2, Q — * ^9 t ne c = 2 daughter 
CFT is obtained. It is interesting to notice that such a daughter CFT gives rise 
to a vortices plasma of half integer vorticity, that is to a fully frustrated XY 
model, as it will appear in the following. 

Its primary fields content can be expressed in terms of a ^-invariant scalar 
field X(z), given by 

X(z) = ±(QU(z) + Q<- 2 \z)), (6) 

describing the continuous phase sector of the new theory, and a twisted field 

<t>{z) = \(QV{z)-QW{z)), (7) 

which satisfies the twisted boundary conditions <\>[e %lx z) — —<f>(z) Such fields 
coincide with the ones introduced in eq. 10} . 

The whole TM theory decomposes into a tensor product of two CFTs, a 
twisted invariant one with c = | and the remaining c = ^ one realized by a Ma- 
jorana fermion in the twisted sector. In the c = | subtheory the primary fields 
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are composite vertex operators V (z) = (z) ip (z) or V q h (z) — U x (z) cr (z), 
where (z) = -j= : e laiX ^ : is the vertex of the charged sector with 

ctl = |,Z = 1, ...,4 for the SU(2) Cooper pairing symmetry used here. 

Regarding the other component, the highest weight state in the neutral 
sector can be classified by the two chiral operators: 

i/>(z) = ^ (: e iV5 *« : + : e <V5 *<-*> : ) , 

$(z) = ^= (: e^*W . _ . e <^(-«) .) . ( 8 ) 

which correspond to two c = | Majorana fermions with Ramond (invariant 
under the Z 2 twist) or Neveu-Schwartz (Z% twisted) boundary conditions |3] in 
a fermionized version of the theory. Let us point out that the energy-momentum 
tensor of the Ramond part of the neutral sector develops a cosine term: 

T^(z) = -i(a0) 2 - I ^cos(2V2^) , (9) 

a clear signature of a tunneling phenomenon which selects a new stable vac- 
uum, the linear superposition of the two ground states. The Ramond fields are 
the degrees of freedom which survive after the tunneling and the Z 2 (orbifold) 
symmetry, which exchanges the two Ising fermions, is broken. 

So the whole energy-momentum tensor within the c = f subtheory is: 

T = T X (z) +T^{z) = -\ {dXf - 1 (d<pf - cos . (10) 

The correspondence with the Hamiltonian in eq. (£Q) is more evident once we 
observe that the neutral current d<j) appearing above coincides with the term 
d<f)— ? of eq. since the ^-term coming there from the frustration condition, 
here it appears in d(f> as a zero mode, i.e. a classical mode. Besides the fields 
appearing in eq. © there are the a (z) fields, also called the twist fields, which 
appear in the primary fields V q h (z) combined to a vertex with charge f . The 
twist fields have non local properties and decide also for the non trivial properties 
of the vacuum state, which in fact can be twisted or not in our formalism. Such 
a property for the vacuum is more evident for the torus topology, where the 
er-field is described by the conformal block (see Section 4). 

The evidence of a phase transition in ladder systems at c = | has been 
investigated in |17) within a CFT framework. Within this framework the ground 
state wave function is described as a correlator of N 2e Cooper pairs: 

<N 2e a\Y[V^(z l )\0>= ]J (3* - z v )Pf ( — ^ (11) 

i=l i<j' = l v ' 

where Pf (^ z .l z ■ ^ = A y Z1 1 Z2 Z3 _ Z4 ■ • • 1 * s the antisymmetrized product over 
pairs of Cooper pairs, so reproducing well known results |18) . In a similar way 
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+ 



(a) 



(b) 



Figure 2: boundary conditions: (a) untwisted; (b) twisted. 



we also are able to evaluate correlators of AT 2e Cooper pairs in the presence 
of (quasi-hole) excitations with non Abelian statistics 19 j. It is now 

interesting to notice that the charged contribution appearing in the correlator of 
N e electrons is just: < N e a\ Y[f=l C/^ /2 (zi)|0 >= l\^i> =1 (zi- z^) 1 / 4 , giving rise 
to a vortices plasma with i? e // = — |5Zi<j=i ^ n \ z i ~~ z j\ a t the corresponding 
"temperature" = — 2, that is it describes vortices with vorticity v = |! 

On closed annulus geometries, as it is the discretized analogue of a torus, we 
must properly account for boundary conditions at the ends of the finite lattice 
since they determine in the continuum the pertinent conformal blocks yielding 
the statistics of quasiparticles as well as the ground state degeneracy. We have 
two possible boundary conditions which correspond to two different ways to 
close the double lattice, i.e. tpff — + tp^' or </?^ +1 ' ) — + Vf , = 1,2. It is not 
difficult to work out that, for the ladder case, the twisted boundary conditions 
can be implemented only on the odd lattice. 

Indeed for / = A the ladder is invariant under the shift of two sites, so 
that there are two topologically inequivalent boundary conditions for even or 
odd number of sites. In the even case the end sites are of the same kind of 
the starting one, while in the odd case a ferromagnetic line corresponds to an 
antifcrromagnetic one. The even (odd) case corresponds to untwisted (twisted) 
boundary conditions, as depicted in Fig. 2. The odd case selects out two 
degenerate ground states which belong to a non trivial topological sector, indeed 
the system may develop topological order in the twisted sector of our theory. 
Let us notice that in the discrete case not all the vacua in the different sectors 
of our theory are connected. For instance the states in the untwisted sector, 
which correspond to the even ladder, are physically disconnected from those of 
the twisted one, which correspond to the odd ladder. 



4 Symmetry properties of the TM conformal 
blocks 

In Section 3 we identified our chiral fields with the continuum limit of 
the Josephson phase defined on the two legs of the ladder respectively 
and considered non trivial boundary conditions at its ends, so constructing a 
version in the continuum of the discrete system. By using standard conformal 
field theory techniques it is now possible to generate the torus topology, starting 
from the edge theory, just defined in the previous Section. That is realized by 
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evaluating the N— vertices correlator 

(n\V a ( Zl )...V a (z N )e 2 ™ TLo \n) (12) 

where V a (zi) is the generic primary field of Section 3 representing the excitation 
at Zi, Lq is the Virasoro generator for dilatations and r the proper time. The 
neutrality condition a = must be satisfied and the sum over the complete 
set of states \n) is indicating that a trace must be taken. For the present paper 
it is not necessary to go through such a calculation but it is very illuminating 
for the non expert reader to pictorially represent the above operation in terms 
of an edge state (that is a primary state defined at a given r) which propagates 
interacting with external fields at z\...Zn and finally getting back to itself. In 
such a way a 2D surface is generated with the torus topology. From such a 
picture it is evident then how the degeneracy of the non perturbative ground 
state is closely related to the number of primary states. Furthermore, since in 
this paper we are interested in the understanding of the topological properties 
of the system, we can consider only the center of mass contribution in the above 
correlator, so neglecting its short distances properties. To such an extent the 
one-point functions are extensively reported in the following. 

On the torus 0] the TM primary fields are described in terms of the con- 
formal blocks of the Z 2 -invariant c = 3/2 subtheory and of the non invariant 
c = 1/2 Ising model, so reflecting the decomposition on the plane outlined in 
the previous Section. The following characters 



Xo( w n\r) 



Xa(w„|t) 



XJgK^r) 

express the primary fields content of the Ising model with Neveu-Schwartz {Z<i 
twisted) boundary conditions, while 

Xo(w n \T)K Q (w c \T) + XI (w n \T)K 2 (w c \T) (13) 

XiK|r) {Kx{w c \t) + K 3 (w c \t)) (14) 

16 

Xk(w n \r)K (w c \T) + Xo{w n \T)K 2 (w c \T) (15) 

represent those of the ^-invariant c = 3/2 CFT. They are given in terms of a 
"charged" K a {w c \r) contribution, (see definition given below) and a "neutral" 
one X@( w n\ T )) (the conformal blocks of the Ising Model), where w c — lnz c 
and w n = hi z n ar e the torus variables of "charged" and "neutral" compo- 
nents respectively. 



1 ( / flaKlr) l 04(w n \T) \ 

2 VV Vir) + V V(r) ) 

l( / g 3 (w»|r) l 9 4 (w n \T) \ 
2 V V V(j) V ^ / 

l e 2 (w n \T) 
V 2 ^ T ) 



c=3/2/ I I \ 

X ( o) {w n \w c \T) = 

c=3/2/ I I \ 

(W n \w c \T) = 

X C { 2) ,2 ( W n\Wc\T) = 
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In order to understand the physical significance of the c = 2 conformal blocks 
in terms of the charged low energy excitations of the system, let us evidence their 
electric charge (magnetic flux contents in the dual theory, which is obtained by 
exchanging the compactification radius R? e — > R^ in the charged sector of the 
CFT). In order to do so let us consider the "charged" sector conformal blocks 
appearing in eqs. (|13I15|) : 



K 2 i+z(w\t) 



1 



r)(r) 



e 



2l+i 
4 





(2w\4r), V(U)G(0,1) 2 



(16) 



They correspond to primary fields with conformal dimensions h 2 i+i 



in- 2 
2 "(/,») 



I (~L±i 



'^ti + 2<5( (+i ) ) 2 and electric charges 2e ^ J (magnetic charges in the dual 

theory || (ay^Rx)), Rx = 1 being the compactification radius. More explic- 
itly the electric charges (magnetic charges in the dual theory) are the following: 



; = 1 = q e = Ae (q m = 2§§ ) 

1 = 1 t = q e = 2e (<? m =|f) 

1 = 1 = 1 q e = e (q m = \%) 

1 = 1 i = l qe = Ze (g m = ftf) 



(17) 



If we now turn to the whole c = 2 theory, the characters of the twisted sector 
are given by: 



X( 0) (w n \w c \T) 

x^C^KIt) 

Xl Q) {Wn\w c \T) 



X±{Wn\T){xl 3/2 (w n \w c \T) +X2 3/2 (Wn\w c \T)j (18) 

X± (xo + Xf) (Ko+K 2 ) 
Xo(wn\r) +Xl(w n \r)) X C i~ 3/2 (wn\w c \T) 



Xi (xo + Xi) (K!+K 3 ) 

X^(w„\t) {xl~ 3,2 (Wn\w c \T) - X2~ 3,2 {Wn\w c \l 

(xo - Xi) (Ka - K 2 ) 

Xo(w n \r) - £1 (m>„|t)) Xl~ 3/2 (w„|w c |r) 



(19) 



(20) 



(21) 



Xi (xo-X|) (K1 + K3). 
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Furthermore the characters of the untwisted sector are [3] : 

X(o)OnK|r) 
X(i)(WnK|r) 



X (0) (w„|w c |r) 



X (1) (w„|w c |r) 



X(o)(w„|w c |r) 

Such a factorization is a consequence of the parity selection rule (m-ality) , which 
gives a gluing condition for the "charged" and "neutral" excitations. The con- 
formal blocks given above represent the collective states of highly correlated 
vortices, which appear to be incompressible. In order to show the corresponding 
symmetry properties it is useful to give a pictorial description of the conformal 
blocks appearing in eq. (|16|) . To such an extent let us imagine to cut the torus 
along the A-cycle. The different primary fields then can be seen as excitations 
which propagate along the -B-cycle and interact with the external Cooper pair 
at point w = ^ hi£. We can now test the symmetry properties of the charac- 
ters of the theory (given above) by simply evaluating the Bohm-Aharonov phase 
they pick up while a Cooper pair is taken along the closed A-cycle. In order 
to do that, it is important to notice that the transport of the "Cooper pair" 
from the upper (with isospin up) leg to the down (with isospin down) leg can 
be realized by a translation of the variables w c and w n , which must be identical 
for the "charged" and the "neutral" sectors. So, under a 27r-rotation the torus 
variables transforms as Aw c — Aw n = 1- It is easy to check that: 

K j(w c + 1\t) = K 0i 2(w c \t), K lt3 {w c + 1\T) = -K lt3 (w c \r). (27) 

The change in sign given in eq. I|27[) under a 27r-rotation is strictly related 
to the presence in the spectrum of excitations carrying fractionalized charge 
quanta. Now, turning on also the neutral sector contribution in the Cooper pair 
transport along the A-cycle we obtain in a straightforward way: 

Xo,i( w « + 1|t) = x ,i(w„|r), x±,{w n + 1\t) = ix^{w n \T) (28) 

the same is true for the characters X/3- A phase i = appears in the Xi 
character due to the presence of a half-flux. 
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Xo(wnk)X(cT) 3/5 \wn\wc\r) + Xi(w„|r)X(^ 3/2 '(w n \w c \T) (22) 
(XoXo + XiXi) K + (xoXi + X|Xo) K 2 
Xo{w„\t)x c { 2) /2 (w„\w c \t) + xi(w n \T)x c { o } 3/2 (w n \w c \T) (23) 



(xoXi + XiXo) K + (xoXo + XiXi) 



K 2 



Xo(w„|t)X( ) (w n \w c \T) - xi (w n \T)x c {2) 3/2 (w n \w c \T) (24) 



K-2 



(XoXo - X|X|) K + (xoXi ~ XiXo) 
Xo(w„|t)x(2) 3/2 (w„K|t) - Xi(wn|r)x^ 3/2 (w„|w c |r) (25) 
(XoXi ~ XiXo) K + (xoXo ~ XiXi) K 2 

X±{Wn\T)x C ( i) /2 {Wn\w c \T) = X^Xi (#1 + K 3 ) . (26) 



As a result the ground state described by eq. (|26[) : 



X(0){wu\w c \t) = X^XJjj {K x + K 3 ) 



(29) 



does not change sign under, Aw c — Aw n — 1, the transport of a Cooper pair 
along the closed A-cycle. In fact the negative sign coming from the continuous 
phase sector being compensated by the negative sign coming from the other 
sector! Of course the same is true for all the other characters of the untwisted 
sector, i.e. we cannot trap a half flux quantum in the hole in the untwisted 
sector. 

Instead in the twisted sector the ground state wave- functions show a non 
trivial behavior. In fact under Aw c = Aw n — 1 



X (0) (w n + M w c+Mr) = +iX( S)) {w n \w c \T), X(i){w n +l\w c +V\T) = -ix^(w n \w c \T), 



The change in phase given above evidences the presence of a half flux quantum 
in the hole. In fact in the twisted case geometry (see Fig2) the Cooper pair 
flows along the ladder and changes isospin in a 27r-period, so implying that in 
such a case the transport of a Cooper pair from a given point w on the ^4-cycle 
to the same point has a 47r-period, that is it corresponds to Aw c = Aw n = 2. 
Under this transformation the characters get the following non trivial phase: 



so explicitly evidencing the trapping of in the hole. 

It is worthwhile to notice that the properties just discussed are indepen- 
dent of the short distance properties of the vortices plasma, the only crucial 
requirement for its stability being the neutrality condition. 

5 Topological order and "protected" qubits 

The aim of this Section is to fully exploit the issue of topological order in a 
quantum JJL. In order to meet such a request let us use the results of the 2- 
rcduction technique for the torus topology 0|. For closed geometries the JJL 
with Mobius boundary conditions gives rise to a line defect in the bulk. So it 
becomes mandatory to use a folding procedure to map the problem with a defect 
line into a boundary one, where the defect line appears as a boundary state. 
The TM boundary states have been constructed in Ref. [5] together with the 
corresponding chiral partition functions and briefly recalled in the Appendix. 
In particular we get an "untwisted" sector and a "twisted" one, corresponding 
to periodic and Mobius boundary conditions respectively. That gives rise to 
an essential difference in the low energy spectrum for the system under study: 
in the first case only two-spinon excitations are possible (i.e. SU (2) 2 integer 
spin representations) while, in the last case, the presence of a topological defect 
provides a clear evidence of single-spinon excitations (i.e. SU(2) 2 half-integer 
spin representations). All that takes place in close analogy with spin-1/2 closed 



(30) 



xf 0A) (w n + 2\w c + 2\T) 



xf 0A) (w n \w c \T), 



(31) 
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(a) 



Figure 3: degenerate ground states for: (a) odd ladder, V — ±1; (b) even ladder, 



zigzag ladders, which are expected to map to fully frustrated Josephson ladders 
in the extremely quantum limit. In this way all previous results obtained by 
means of exact diagonalization techniques |2()j are recovered. 

In the following we discuss topological order referring to the characters which 
in turn are related to the different boundary states present in the system through 
such chiral partition functions p^. We can build a topological invariant V = 
l\ Xp where 7 is a closed contour that goes around the hole. Such a choice allows 



us to define two degenerate ground states with V = 1 and V = — 1 respectively, 
labelled |0) , |1), as in Fig. 3(a): there must be an odd number of plaquettes (odd 
ladder) to satisfy this condition which is imposed by the request of topological 
protection. The value V = — 1 selects twisted boundary conditions at the ends 
of the chain while V = 1 selects periodic boundary conditions (see Fig. 3(b)). 
The case of finite discrete systems has been discussed in detail in Ref . [21] , our 
theory is its counterpart in the continuum. 

In fact it is now possible to identify the two degenerate ground states shown 
in Fig. 3(a) with the twisted characters I|18I19|I : 



where it is meant that w n is fixed to zero (see also (I55[l and <|56l) of Appendix), 
then it is possible to remove such a degeneracy through vortices tunneling. The 
last operation is also needed in order to prepare the qubit in a definite state. 

In fact such two states can be distinguished from all the other ones. First of 
all, being in the twisted sector they trap a half flux quantum inside the central 
hole (see 1)3111). The trapped half flux quantum can be experimentally detected, 
so giving a way to read out the twisted states of the system. 

Now, upon performing an adiabatic change of local magnetic fields which 
drags one half vortex across the system, i.e. through the transport of a half flux 
quantum around the £?-cycle of the torus, we are able to fully distinguish such 
two states with respect to the other pair of states, x^o) an d 1 °f the twisted 



V = l. 



1 



|o)^xJ)Mt). |i)^X(i)K|r), 



(32) 



sector. 



In fact, the transformations 



xfo)(^/2kc + T/2|r) 
x5)(7-/2|«,c+t/2|t) 



-ie-^q-^x+^wclr) 
ie-^q-^x+^wclr), 



(33) 
(34) 
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X m (r/2\w c + T/2\T) = 
X 7 1) {t/2\wc + t/2\t) = 



(35) 
(36) 



just imply a flip among the two ground states of the system 



|0)-|1) and |1)-|0), 



(37) 



while the other two twisted states are annihilated. 

Finally it should be noticed that the presence of the twist operators a, de- 
scribed here by the Ising character gives rise to non Abelian statistics, 
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which can be evidenced by their fusion rules a a = 1 + ip |18|[4*]. 

6 Brief summary with comments 

In this contribution we presented a simple collective description of a ladder of 
Josephson junctions with a macroscopic half flux quanta trapped in the hole of a 
ladder. It was shown how the phenomenon of flux fractionalization takes place 
within the context of a 2D conformal field theory with a Z2 twist, the TM. 
The presence of a Z2 symmetry indeed accounts for more general boundary 
conditions for the fields describing the Cooper pairs propagating on the ladder 
legs, which arise from the presence of a magnetic impurity strongly coupled with 
the Josephson phases. It should be noticed that it would be useful to extend our 
approach to a generic frustration / = ^ • For closed geometries and in the limit 
of the continuum the phase fields tp( a ^ defined on the two legs were identified 
with the two chiral Fubini fields of our TM, and a correspondence between 
the energy momentum density tensor for such fields (or better the X and <j> fields 
of eqs. ©-0) and the Hamiltonian of eq. traced. For such geometries it 
was also indicated that the Kosterlitz-Thouless vortices were recovered. 

Furthermore it was shown that for closed geometries the JJL with an im- 
purity gives rise to a line defect, which can be turned into a boundary state 
after employing the folding procedure. That enabled us to derive the low en- 
ergy charged excitations of the system as provided by our description, with the 
superconducting phase characterized by condensation of 4e charges and gapped 
2e excitations. Finally, by simply evaluating a Bohm-Aharonov phase, it has 
been evidenced that non trivial symmetry properties for the conformal blocks 
emerge due to the presence in the spectrum of fractionalized flux quanta ^ (ff ) ■ 
As it has been explained before, that signals the presence of a topological defect 
in the twisted sector of the TM. The question of an emerging topological order 
in the ground state together with the possibility of providing protected states 
for the implementation of a solid state qubit has been also addressed pQ . Notice 
also the different behavior under transport of the 2e and 4e excitations as it is 
well evidenced by the Bohm-Aharonov phase. Indeed while the transport of a 
2e excitation along the A-cycle induces a —1 phase factor, in the 4e excitation 
transport the phase factor is trivial That is a consequence of the symmetry 
of the 4e field with respect to the leg index. 
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Finally we have shown that Josephson junctions ladders with non trivial ge- 
ometry may develop topological order allowing for the implementation of "pro- 
tected" qubits, a first step toward the realization of an ideal solid state quantum 
computer. Josephson junctions ladders with annular geometry have been fab- 
ricated within the trilayer Nb/Al — AlO x /Nb technology and experimentally 
investigated |22| . So in principle it could be simple to conceive an experimental 
setup in order to test our predictions. 

It is interesting to notice that the presence of a topological defect has been 
experimentally evidenced very recently for a two layers quantum Hall system, 
by measuring the conduction properties between the two edges of the system 

m 

Appendix 

Here we recall briefly the TM boundary states (BS) recently constructed in [S]. 
For closed geometries, that is for the torus, the JJL with an impurity gives rise 
to a line defect in the bulk. In a theory with a defect line the interaction with 
the impurity gives rise to the following non trivial boundary conditions for the 
fields: 

v4 Q) (x = 0) = T^ Q) (x = 0) - tpo, a = 1, 2. (38) 

In order to describe it we resort to the folding procedure. Such a procedure is 
used in the literature to map a problem with a defect line (as a bulk property) 
into a boundary one, where the defect line appears as a boundary state of 
a theory which is not anymore chiral and its fields are defined in a reduced 
region which is one half of the original one. Our approach, the TM, is a chiral 
description of that, where the chiral cj> field defined in (-L/2, L/2) describes 
both the left moving component and the right moving one defined in (-L/2, 
0), (0, L/2) respectively, in the folded description jS]- Furthermore to make a 
connection with the TM we consider more general gluing conditions: 
4> L {x = 0) = t4>r{x = 0) - ip 

the — (+) sign staying for the twisted (untwisted) sector. We are then allowed 
to use the boundary states given in (J.O for the c = 1 orbifold at the Ising 2 
radius. The X field, which is even under the folding procedure, does not suffer 
any change in boundary conditions Let us now write each phase field as 
the sum ip^ (x) = ip^ (x) + (p^' (x) of left and right moving fields defined 
on the half-line because of the defect located in x = 0. Then let us define 
for each leg the two chiral fields ipi a l (x) = ip^ (x) ± (p^/ (—x), each defined 

on the whole x— axis jJS]. In such a framework the dual fields (p^ (x) are 
fully decoupled because the corresponding boundary interaction term in the 
Hamiltonian does not involve them they are involved in the definition of 

the conjugate momenta H( a ) — (^dxip^^j = y - 3 (a) ^ present in the quantum 

Hamiltonian. Performing the change of variables tpe — X + 0, p)f^ = X — cp 
((p^ — X + (f>, tp^ = X — 4> for the dual ones) we get the quantum Hamiltonian 
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but, now, all the fields are chiral ones. 

It is interesting to notice that the condition 1|38|) is naturally satisfied by the 
twisted field <j) (z) of our twisted model (TM) (see eq. Q). 

The most convenient representation of such BS is the one in which they 
appear as a product of Ising and c = | BS. These last ones are given in terms of 
the BS | a > for the charged boson and the Ising ones |/ >, | f>, | J.> according 
to (see ref. fZT\ for details): 

|X(of 72 >= |0 > ®| T> +|2 > ®| |> (39) 

l41 3/2 >=^74(l 1 >+l 3 >)®l/> ( 4 °) 
| X ( 2 f 72 >= |0 > ®| |> + |2 > ®| |> • (41) 

Such a factorization naturally arises already for the TM characters 0]. 

The vacuum state for the TM model corresponds to the X(o) character which 
is the product of the vacuum state for the c = | subtheory and that of the Ising 
one. From eqs. (|22I24|I we can see that the lowest energy state appears in two 
characters, so a linear combination of them must be taken in order to define a 
unique vacuum state. The correct BS in the untwisted sector are: 

IX((o,o),o) >= -j= (lx+o) > +lxf ) >) = V^(|0 > ®| T f > +|2 > ®| 4 f >) 

(42) 



X((Q,o),i) >= -j= > -1% >) = V2(|0 > > +|2 > ®| T I >) 



(43) 



X((i,o),o) >= 7/| (|X(i) > +lx w >) - V2(|0 > ®|i f > +|2 > ®| T t >) 



(44) 



IX((i,o),i) >= > -|xf X ) >) = \/2(|0 > ®| T I > +|2 > ®||I >) 

(45) 

IX(o)(^o) >= (|1 > +|3 >) ® \D {<po) > (46) 

where we also added the states |x(o) (fo) > m which \Do(<po) > is the continuous 
orbifold Dirichlet boundary state defined in ref. For the special ip = n/2 

value one obtains: 

l*(0) ^4 (Il>+|3 >)«>!//>■ (47) 
For the twisted sector we have: 

| X + > = (|0 > +|2 >) ® (| T / > +|i / >) (48) 

Ixfi) > = (|1 > +|3 >) ® (|/f > +1/1 >)■ (49) 
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Now, by using as reference state \A > the vacuum state given in eq. (|42[l. 
we compute the chiral partition functions Zab where \B > are all the BS just 
obtained [S]: 



<X((o,o),o)IIX((o,o),o)> = X((0,0),0) (50) 

<X((o,o),o)IIX((i,o),o)> = X((1,0),0) (51) 

<X((0,0),0)IIX((0,0),i)> = X((o,o),i) (52) 

<X((o,o),o)IIX((i,o),i)> = X((i,o),i) (53) 

2 '<X«o,o),o)IIX(o)> = X(0) (54) 

Z <X( ( o,o),o)llx+ 0) > = X U ( 55 ) 

z <x((o,o ) ,o)iixi l 1) > = x tr ^ 



So we can discuss topological order referring to the characters with the im- 
plicit relation to the different boundary states present in the system. Further- 
more these BS can be associated to different kinds of linear defects, which are 
compatible with conformal invariance 5 . 
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